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Collective excitations in graphene monolayer are studied. Equations describing collective properties of 
electrons in graphene are obtained. The basic ideas of the method of many-particle quantum hydrodynam- 
ics are used for the derivation. As starting point of the derivation we use the Dirac equation for massless 
electrons which is usually used for description of electrons in graphene [D. E. Sheehy and J. Schmalian, 
Phys. Rev. Lett. 99, 226803 (2007)], where the Coulomb interaction is taken into account. We study 
dispersion properties of collective excitations by means derived here graphene quantum hydrodynamics 
equations (GQHD). We consider graphene in the external magnetic field which directed at an angle to the 
graphene sample. We do it in a linear approximation of the GQHD equations. We observe that the mag- 
netic field directed perpendicular to the graphene plane had no influence on dispersion of the collective 
excitations. For the magnetic field directed at an angle to the graphene we obtain dependence of wave 
dispersion on system parameters: strength of magnetic field, wave vector, direction of wave propagation 
relatively to the magnetic field. 



I. INTRODUCTION 

Unusual properties of graphene conductivity |U], 
[Q] have led to its widespread use in semiconductor 
geterostructures. Important characteristic of three dimen- 
sional (3D) and low dimensional conductors and semicon- 
ductors is a dispersion of collective excitations, particularly 
the plasma waves or plasmons. Knowledge of the plas- 
mon dispersion allows us to analyze various processes in 
geterostructures. A model describing the collective prop- 
erties of graphene is constructed in this paper. We use this 
model for definition of dispersion of the quantum collective 
excitations in graphene at the presence of the external mag- 
netic field. The background of this research is the model 
presented in Ref. yfl. This model describes the dynamic of 
electrons in graphene and reflects essential characteristics 
of the conduction electrons in graphene. The basic equa- 
tion describing microscopic dynamics of the conduction 
electrons in graphene is the many-particle equation which 
can be presented in the form analogous to the Schrodinger 
equation (or the Dirac equation) y|] 



ihd t ip = Hip, 



(1) 



where H is the Hamilton operator. In the absence of ex- 
ternal fields and interaction between particles the Hamilto- 

nian has form H = J2iLi v F^tP'i^ nere N is the num- 
ber of particles in the system, Vp is the Fermi velocity 
of electrons in graphene, pf is the momentum operator, 
pf = —ihdf, df is the derivative on coordinate of i- 
th particle, af are the Pauli matrixes. Described model 
corresponds to the massless behavior of the conductivity 
electrons in graphene 0]. The Coulomb interaction be- 
tween the electrons in graphene is considered in Ref. JH. 
Equation £Q as the many-particle Schrodinger equation is 
not always suitable for description of the collective prop- 
erties in many-particle systems. This problem connects 
with the fact that equation (Q} determines the wave func- 
tion in the 3N (2N, for two dimensional system) dimen- 



sional configuration space, whereas the collective process 
realize in the 3D (2D) physical space. Therefore, it is im- 
portant to construct a method in the 3D (2D) physical space 
H. The method of quantum hydrodynamics (QHD) solves 
the problem of transition from configuration space to phys- 
ical space. This method was suggested for many-particle 
system and has been developed for a wide class of physi- 
cal systems. Method of the many-particle QHD (MPQHD) 
allows us derive equations for quantum observable values 
evolution in the 3D (2D) physical space. Equations of con- 
tinuity, momentum balance, energy balance, momentum 
evolution (for particles with spin) and polarization (for po- 
larized particles) are appeared from the Schrodinger equa- 
tion. Most of them are analogous to classic hydrodynamic 
equations. In connection with described above the method 
is called "Quantum hydrodynamics". Method of the QHD 
is a modern powerful method for studyin g of the collective 
properties in systems of charged llql- iflOTI -. and neutral 
ill ill . lfl3ll particles (see also review lll4ll ). Unusual prop- 
erties of electrons of conductivity in graphene mapped in 
the equation CO and mapped in the form of the QHD equa- 
tions for graphene (GQHD), which we present below. The 
equations of the GQHD in this paper has the extreme dif- 
ference from the QHD equations obtained previously |@] 
- lfl4ll . In context of using of the Dirac equation for de- 
scription of graphene electrons we note that in Ref. lfl5ll 
the QHD equation were derived from the bispinor Dirac 
equation for massive particles. Obtained there equations 
were averaged on ensemble for receiving of equations for 
many-particle system. Comparison of the QHD equations 
for graphene with the usual QHD equations ifioll . [ 14 1 we 
present below during the derivation of the GQHD equa- 
tions for graphene. As expected unusual properties of elec- 
trons of conductivity in graphene lead to exotic spectrum of 
elementary excitations. In 2D electron gas (2DEG) the ba- 
sic type of the collective excitations is the Langmuir waves, 
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the dispersion dependence of these waves are 

2ire 2 n k 
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where e, m are the charge and mass of electrons (e we 
consider as algebraic quantity, e = — |e|, for electrons), 
n a is the 2D concentration of electrons, k is the absolute 
value of the wave vector. For the graphene electrons in the 
absence of external fields the dispersion equation is 
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So, we did not obtain the contribution of interaction in 
dispersion, and, consequently, from our description the 
graphene Langmuir frequency (see Ref. [1] p. 414) does 
not follow. 

The recent studies of excitations in graphene being in a 
magnetic field have employed the picture of single electron 
excitation and it's transition between energy levels, such as 
the Landau levels, which exist for single charged particle 
in an external field. Consideration of microscopic picture 
of excitations show that some of these excitations lead to 
the formation of collective excitations 

If we want to study macroscopic collective excitations 
in the graphene, we do not need to consider one electron 
excitation. For described purpose we can use hydrody- 
namics or kinetics methods as it has been used for two 
dimensional electron gas, usual three dimensional carri- 
ers and semiconductors. These method get used to be 
used for getting of characteristics of the collective excita- 
tions in described objects, such as the Langmuir waves, the 
magneto-sound waves, the Bernstein modes, and etc. So, 
we can see that obtaining of the hydrodynamics descrip- 
tion of graphene is the important task. The method of the 
many-particle QHD gives us possibility to derive GQHD 
equations, which make more rich a set of theoretical tools 
of graphene excitations studying. 

Our paper is organized as follows. In Sec. II we present 
derivation of the GQHD. In Sec. Ill we describe the 
method for calculation of the dispersion of the collective 
excitations. In Sec. IV dispersion of collective excitations 
for the graphene electrons in the external magnetic field is 
studied. In Sec. V the brief description of obtained results 
is presented. 



II. CONSTRUCTION OF THE MODEL 

We derive basic equations by means of the MPQHD 
method @], ifToh . We use the many-particle spinor mass- 
less Dirac equation (3], 0] 

\ i ' i,j=£i / 

(2) 

The following designations are used in the Hamilto- 
nian ©: Df = -iHd? - e^Jc, ip hext , A" 



- is the potentials of the external electromagnetic field, 
E t ,ext = -V<Pi,ext ~ d t K,ext is the electric field, B hext = 
curlAi >ex t is the magnetic field, quantities e i; m^are the 
charge and mass of particles, h-is the Planck constant, and 
Gij = - is the Green functions of the Coulomb in- 

teraction. Replacing —ihd" by Df were used in Ref. rf20ll 
(see. p. 127). It was made for account of external magnetic 
field In equation (0 the spinor wave function ip = ip{R, t) 
depend on 2N coordinates R = [r 1,..., r N ] and time, 
where r< = [», , y$] is the 2D coordinates of each particle. 
Potentials <p i>ext = Lp ext (r t ,t), A<* ext = A^r^t) also 
depend on 2D variables. This fact has deep consequences. 
Potential part of electric field connected with the scalar po- 
tential via space derivative: Ej = — Vj^. Consequently 
in equation (0 there is no contribution of external electric 
field directed perpendicular to the graphene plane E z (Con- 
tribution of E z might appear via d t A z ). Physically, there 
is no limitation on attendance of z projection of electric 
field and it's action on graphene electrons. Especially if the 
graphene sample is the_part of the geterostructure or spin- 
field-effect transistor 111, UM, Q, there exist the 



contribution of external electric field in normal direction to 
the graphene plane. The magnetic field vector to be 

B = curlA = e x (d y A z — d z A y ) 

~\~£y{d z A x — d x A z ) + e z (d x A y — d y A x ), 

two component of the vector potential of the magnetic field 
A x , A y are presented in Hamiltonian (O, and they does 
not depend on coordinate z. Therefore equation (O con- 
tain z component of the magnetic field only. In this paper 
we interested in action of the external magnetic field di- 
rected at angle of graphene plane. Therefore, we general- 
ized GQHD equations including whole vector of magnetic 
field B ■ e z ^ B = [B x ,B y ,B z ). 

We assumed that the Pauli matrices satisfy the following 
commutation relation: 



(3) 



Graphene is the 2D structure and electrons of graphene 
are located in the plane. As we describe above, in 2D case 
the electrons has two coordinate x and y, but spin of the 
electrons can be directed in all direction, particularly, in di- 
rection of z axes, perpendicular to the graphene plane. This 
fact is accounted by formula (O. Two projection of spin 
operator are contained in the Hamiltonian these are a x 
and (Ty. Spin operator projection on z axis is appeared dur- 
ing derivation of the GQHD equations due to commutation 
relation (O. 

Now, we present the derivation of model of the electron 
collective dynamics in graphene. As we mention above 
we describe the derivation of the MPQHD equations from 
equation (O. The first step is we present the definition for 
density of probability for the conduction electron system in 
the physical space. 

n(r,t) = J2 J dRj2^r-r t )r(R,t)^(R,t) (4) 



3 



where dR = Ylp=i dr p - 

The quantity n(r,t) can be considered as the 2D con- 
centration of the conductivity electrons. For studying of 
the time evolution of the concentration we differentiate the 
concentration (0 with respect to time and use equation (0. 
In the result, we receive an equation which has the form of 
the continuity equation: 

d t n{r,t) + Vj(r,t) = 0, (5) 

where 



x UC(J2,i) o-fV (R,t) + h.c. 



(6) 



and j(r, t) = v F S(r, t). 

The quantity S(r, t) describes the spin density of the par- 
ticle system. Consequently, we have equation 

d t n(r,t) +v F VS(r,t) = 0. (7) 

Coordinate vector r has only two component. Conse- 
quently, equation (0 contains two component of the spin 
density vector S, these are S x and S y . Our next step in con- 
struction of the model of collective motion is obtaining of 
an equation for spin evolution. For this aim we differentiate 
the quantity S(r, t) with respect to time and use equation 
0. Because we known the third component of spin den- 
sity vector we can study the evolution of all components 
of this vector. Therefore, we derive equation of evolution 
for S = [S x ,S y , S z ]. On this way we have equation of the 
spin evolution: 

d t S a (r,t) + v F d a n(r,t) = ~s a ^J^(r,t) (8) 



Here new physical quantity is arisen: J^f (r, t). The evi 
dent form of J"f (r, t) is 



dRj2$(r 



,1 
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the usual QHD 0]- OH. In the usual QHD as in the clas- 
sic hydrodynamics the first equation is the continuity equa- 
tion where described the changing of number of particles in 
vicinity of the point of the physical space in consequence 
of the particles current. Instead of that we obtain the con- 
nection of particles number changing and the spin density 
0. 

We have obtained the equation (0 instead of the mo- 
mentum balance equation (the Euler equation) in the usual 
QHD. The last one accounts the influence of particles in- 
teraction as each other and with the external fields. Equa- 
tion dSJ does not contain information about the interaction. 
Equation (0 differs from the usual equation of the spin 
evolution (the Bloch equation) which contains the vector 
product of the spin density and magnetic field whose leads 
to the spin evolution. 

At this step we have no closed set of two equations (0 
and (0. Thereto, received equations do not contain in- 
formation about interaction. Thus, basic model, presented 
by equation (0, contains the Coulomb interaction between 
electrons. Therefore, we have the series of equations which 
is an analogous to the BBGKI series lE^ll . [25] in the clas- 
sic physical kinetics. We can obtain next equation of the se- 
ries, namely equation of evolution for the tensor «/^f ( r > 
This equation to be 

d t Jtf(r,t)+v F d a J (r,t) 
-hv%e a ^d p d-<S 5 {r,t) = --v 2 F e a ^n{r, t)£ 7 (r, t) 



2v 2 

+— F -e a ^W^(r, t) + ev F S a (r, t)E p (r, t) 
n 

- e 2 v F S a (r,t)8 j dr'G(r,r')n(r',t). (10) 

Equation (fTOl contains interaction and two new quantities, 
these are J Q (r, t) and n a ^ 7 . They are 

J a (r,t) = v F Y / j dRj2S(r- ri )^x 



(R,t) + h.c. 



(9) 



This is the tensor of spin current. We study the 2D struc- 
tures existing in 3D space. It leads to existence of z com- 
ponent of vectors for several physical quantities, as for 
the spin density vector. The spin current J^f(r,t) is de- 
fined via two-vectors (|9]). These are the vector of spin a a 
and derivative vector operator , the last one is the vec- 
tor with two components. Therefore, we have the tensor 
J^f (r, t) whose first index take three value x, y, z, but the 
second index take two value x, y. We have obtained two 
equation of the GQHD, it is equations (0 and (0. These 
equations significantly vary from the first two equations of 



and 



x \r a {RM D ^) (R,t) + h. c . (id 



x U;(R, t) lofDfDW) (R, t) + h.c.J . (12) 

Equation (fTOl is the first equation in system of the 
GQHD equations which contain the magnetic field. In 
equation ( fTOl) we have made generalization and consider all 
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three component of the magnetic field instead of B z which 
presented in the Hamiltonian (0. 

The second rank tensor e al3 ' y B'J xt contains strength of 
the external magnetic field and gives contribution to the 
evolution of spin current. B x and B y do not include in 
equation but they can make influence on particles dy- 
namic. 

Quantity J(r, t) is an analog of the current of probabil- 
ityor the momentum density in the usual QHD |@], lfl4fl . 
II lfjll . This is very important physical quantity and we pro- 
pose that our model must be complete by one more equa- 
tion - the momentum balance equation. For derivation of 
the momentum balance equation, i.e. equation of evolution 
of J(r, t), we apply the same method as used above. We 
differentiate J(r, t) with respect to time and use the equa- 
tion (0. In the result we have 

d t J a (r,t) +v F d JP a = ^Le^S^B^ 

c 

+ ev F nE a - e 2 v F nd a J dr'G(r, r>(r', t). (13) 

Equation [\3[ is very similar to the Euler equation in the 
usual QHD. But instead of the Lorentz force we have vector 
product of the spin density and the magnetic field. 

Equation ( fT3t . just as equation (fTot . contains interaction 
between particles. In fact equation (fl~3T > is an analog of the 
Euler equation in the usual QHD fl, 0], Efl]. Ob- 
tained set of equations 0, (HJ, ([Tot and dTSl is not close 
because the set of equations contains quantity II a/37 (r, t). 
We can obtain equation for iI Q ^ 7 (r, t), but in this way we 
will obtain new physical quantities. We need to make a 
closed set of the equations (0, dU), (fTot and O, express- 
ing iI Q ^ 7 (r, t) via physical quantities included in these 
equations. For closing the set of equations 0, ([8]), ([TUt 
and (fl"3l l we use method of approximate calculations of the 
hydrodynamic variables developed in |@], ifTTll . Following 
the paper |@], Hill we obtain 

n^(r,t) = ^(r,t)^M^M 

n (r, t) 

+ (r,t)-h 2 S a (r,t) g!gygpl . ( 14) 

y/n[r,t) 

where £> Q/37 (r, t) describes the contribution of the thermal 
motion, so, it is an analog of the tensor of kinetic pressure 
in the usual hydrodynamics. In this paper we do not ac- 
count the thermal motion. Consequently we suggest that 
£>"° 7 (r, t) is equal to zero. The last term in formula (fl~4l) is 
an analog of the quantum Bohm potential (see for example 

E3). 

Now, we have a closed system of the equations 0, ([8]), 
(Mob . (PT3T l and £Q). Using the GQHD equations ©, ®, 
( fTUl) . (fT3l and (fT4l we can study the collective excitations 
and they properties, i.e. dispersion dependence and incre- 
ments of instabilities. 
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FIG. 1. The figure presents the picture of the self-consistent in- 
teraction in the system of charged particles. 

Developing method allows us to derive the energy bal- 
ance equation which give us possibility to study the influ- 
ence of the temperature on the graphene dynamics. 

Self-consistent field approximation: discussion 

In last terms of equations ( fTol) and ( fT3l) we have used 
the self-consistent field approximation. Here we present 
general expressions for two-particle functions appeared in 
equations dTOl) and ( fT3l . and also describe the meaning of 
the self-consistent approximation. This approximation was 
suggested by A. A. Vlasov in 1938 12611 for many-particle 
system of charged particles. 

The Fig. [Qpresents the picture of the self-consistent in- 
teraction between charges particles. Total charge of the 
region 2 interacts with the total charge of the region 1. 
Changing the extreme point of a radius vector (or shift- 
ing the region 2) we can scan whole space. In this way 
we obtain action of the external charges on the region 1. 
Changing position of the region 1 and repeating described 
operation we obtain action surrounding charges on each re- 
gion of space. This is a picture of the self-consistent inter- 
action in fixed moment of time and this picture governs an 
evolution of particles in system. This picture of interaction 
is typical for the classic physics, where we need to obtain 
smooth functions describing the collective motion. For that 
is necessary to average at physically infinitesimal volume 
(sketched circle). In the quantum mechanics, the concen- 
tration, the spin density, the current density, etc, are defined 
via wave function and we can consider described picture on 
interaction of separate particle instead of the space regions. 

In general case the last terms in the equations ([TUt and 
(fT3l) contain the two-particle function which in the self- 
consistent field approximation expressed via one particle 
functions and in this approximation we obtain a closed set 
of equations. In the equations ([TUt and (fT3l l the following 
two-particle functions were appeared. Two-particle con- 
centration 

n 2 (r, r',t) = ^2 / ^Rx 
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x Y / 6(T-r i )6(r , -T j )1>+(R,t)il> a (R,t) (15) 

is the average of the product of the concentration operators 
SiLi <H r ~~ r i)' where A?" is the number of particles in the 
system. The quantity 



yields the homogeneous system of algebraic equations 



x J2 S ( r - r iW- r J>^t(R,t)[a^) (R,t) 

(16) 

is the two-particle function of the concentration and the 

spin density X^Li a i^( r ~ r i)- Technically, the self- 
consistent field approximation corresponds to the factor- 
ing of the two-particle functions in product of one-particle 
ones: 



n 2 (r,r',i) ->■ n(r, t)n(r', t) 



(17) 



and 



j 2 a (r,r',t) ^v F S a (r,t)n(r',t). (18) 

It was used at derivation of the equations ( PTOl) and (fT3l) . 
This approximation is suitable at consideration of the long- 
range interaction, particularly for the Coulomb interaction. 

If we interest in more detail description of interaction we 
need to consider many-particle correlations. For example, 
for the two-particle concentration the correlation to be 

g(r, r', t) = n 2 (r, r', t) - n(r, t)n(r', t), 

where g(r, r', t) includes a quantum correlation caused by 
exchange interaction. A method of correlation calculation 
was developed in Ref.s JfJ, ifTlll and l27ll . This method 
also can be used for the graphene description, and for it's 
further and more detailed studying. 



III. METHOD OF CALCULATION OF WAVE 
DISPERSION 

We consider the small perturbation of the equilibrium 
state like 

n = n Q + Sn, S a = S% + 5S a , S || B 

J° = + 5 J a , J A f = + 6J$ , (19) 

where B is the external magnetic field, B = 

[B Ox ,0,B Qz ]. 

Substituting these relations into system of the GQHD 
equations CJ, ©, ( fTOl ), ( fl"3l and (TT^t and neglecting by 
nonlinear terms, we obtain a system of linear equations 
with constant coefficients. 

Passing to the following representation for small pertur- 
bations 5f 

5f = f(u), k)exp(—iujt + ^kr) 



(20) 



iu5S a + iv F k a 5n = -^e aM b f K ] , (21) 



iu5n + iv F k5S = 0, 
2 

T 

,2 



-ico5J a + iv F k?6J 0a = ^e^BZSS* 

c 

+ ik a e 2 v F n (27r/k)5n 

-iujSJI? + iv F k a 5J + hvpk^k^e^dS 5 

= --v 2 F e a ^B^5n + v F he a " u —S^kHn 
c n 

+ ie 2 v F S£kP(2Tr/k)5n. (23) 

We assume that the spin density magnitude has a nonzero 
value. Expressing all quantities entering in the set of equa- 
tions in terms of the spin density, we get the equation 

A Q ^(u;,k) • S p (u,k) = 0, 

where 

A Q/3 (w,k) = uj 2 5 aP -v 2 F k a k? 
+2-—e a ~ ai e ,&v B"k fi k s k i1 —^- 



(22) 



and 



c hu 



uj 2 — vpk 2 



-l—e^^e^BlW "" F , + 2e a ~<»e^»vlk~ < k s 
he oo 2 — vi.k 2 



UJVl 



-2-£ a ^^e^^B v —vlkP+2e ai5 £ hlv v F kPk~<V l S v — 
c hoj n u) 



+ 47ne 2 £ ai5 kPk~«S s 



1 



Hub 



(24) 



Contribution of the Coulomb interaction is presented by the 
last term in the dispersion matrix A a ". 

In the absence of external magnetic field the dispersion 
matrix A Q/3 can be presented as 

A Q/5 (w,k) = uj 2 5 aP - v 2 F k a k p 

+2e ai »e f}S »v 2 F Wk s . 
Dispersion equation to be 

de£A(w,k) = 0. (25) 
This equation can be represent in the form 

(A xx A yy - A xy A yx )A zz = 0, (26) 
and splits into two 

A xx A yy A xy A yx 0, (27) 

and 

A zz = 0. (28) 
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IV. GRAPHENE IN THE MAGNETIC FIELD 

In this section we consider dispersion properties of 
waves in graphene placed in the external magnetic field 
which parallel or perpendicular to the XY plane, we sup- 
pose that graphene is located in the XY plane. 

We begin this chapter with the consideration of equation 



1.04 



^-xy^-yx 



0. 



(29) 



For the magnetic field perpendicular to the graphene 
plane B = B z e z from (l29l we have 

uj = v F k. (30) 

It is the same result as in the absence of the external field. 

Suggesting what Sq x ~ B 0x we can admit that equation 
(l29l contains terms proportional to B 0x and Bq x . Consid- 
ering the limit of approximately small in-plane magnetic 
field B 0x we can neglect by the terms proportional to the 
Bq x . In the result we find an equation which is the equa- 
tion of eight degree of uj. However, at the considering val- 
ues of the parameters of the system we get that this equa- 
tion has two solutions only. These solutions reveal strong 
nonmonotone dependence of the collective excitation fre- 
quency uj on the angle between direction of wave propaga- 
tion and the direction of the in-plane magnetic field (|2j. 

From equation ( |29l we find that in described approxima- 
tion £ depends on three dimensionless parameters 



eB, 



Ox 



chk 2 



<0, 



e 2 S 0z 



X — Sqx/Uq, 

and angle of wave propagation 9. The angle 9 is defined as 
cos 9 = k x /k and sin 9 = k y /k. 

Let's consider equation A zz = 0. Dispersion equation 
has following explicit form 



uj 2 + 2vlk 2 



e uVpkBo 



he uj 2 

or in dimensionless form to be 

.2 , „ 2a£ 



Vpk 2 







? + 2 



where 



and 



UJ 



kv r 



e_kB = _\e[Bo cos6 
he k 3 he k 2 



(31) 



(32) 



(33) 



(34) 



£ describe the frequency of the wave divided by kv F , 
a present the contribution of the external magnetic field 




-3 -2 -1 



FIG. 2. The figure presents the dependence of the reduced fre- 
quency £ (33t , appearing from equation ((29} in presence of the 
in-plane magnetic field Bq x , on the angle 0, which is the angle 
between axes Ox and direction of the wave propagation. Dimen- 
sionless parameters of the system are choose to be equal to the 
following quantities « = — 10 -6 , x — 0.1, and fj, = 10 -3 . Blue 
curve (upper) lying at £ > 1 corresponds to one of solutions of 
equation ( |29l l and red curve (lowest) lying at £ < 1 described the 
second solution. 



which parallel to the graphene plane and 9 is the angle be- 
tween the magnetic field and the direction of the excitation 
propagation. 

In the absence of the external field or if the magnetic field 
directed perpendicular to the graphene plane from d32l l we 
have 

£ 2 + 2 = 0, 

thus the equation A zz = has no wave solution. 

Here we consider equation (l32l l in the case when the 
magnetic field parallel to the graphene plane, then the equa- 
tion (l32t we represent in the form 

£ 4 + C - 2a£ - 2 = 0. (35) 

The dependence of £ on a is presented Fig. [5J This solution 
shows no instabilities. We obtain this solution from the 
equation 

A ZX S X + A zy S y + A ZZ S Z = 0, (36) 

since A zz = and S z ^ we can note that obtained wave 
solution contains wave of spin, where amplitude of spin 
density directed perpendicular to the graphene plane and, 
thus, to the direction of the wave propagation. We need 
to pay attention to the fact that the Coulomb interaction 
contained in the basic equations and the dispersion matrix 
does not influence on found dispersion dependence. 



V. CONCLUSION 

We investigated the influence of the external static uni- 
form magnetic fields on the dispersion properties of lin- 
ear waves of the graphene electrons. We supposed that 
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FIG. 3. The figure presents the dependence of the reduced fre- 
quency £ d3 3 b on parameter a d34b . 

the external field directed at an angle to the plane where 
graphene is located. We paid attention to the particular 
cases when the angle between the magnetic field direction 
and the graphene plane equal to or tt/2. In the absence of 
external fields the dispersion dependence of the collective 
excitation has form u) = Vpk (see formula (l30ll). If the 
magnetic field perpendicular to the graphene plane there 
is no changes in the dispersion in compare with the case 
of graphene in the absence of the external field. We ob- 
tained dispersion dependence for the case magnetic field 
parallel to the graphene plane and studied dependence of 
the frequency of the collective excitations on strength of 
the external magnetic field and angle between the magnetic 
field and the direction of excitation propagation. We also 
showed that in the presence of an in-plane magnetic field 
B 0x , along with B 0z , solution to = Vpk splits on two. 

For studying of described problem we derived system 
of the QHD equations for electrons in graphene. For this 
derivation we used the method of the MPQHD. Obtained 
GQHD equations consist of four equations: equations of 
concentration evolution, spin evolution, current evolution 
and spin current evolution. In this equation we made gen- 
eralization and included in the equations the contribution 
of B x , By along with the B z . 
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